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ABSTRACT 



The separation of precipitation from non-precipitation 
events using a matrix of GOES-S digitized infrared and 
visual satellite data was studied. Precipitation verifica- 
tion was conducted with collocated surface observations. 
The data set consists of 70,623 surface observations, of 
which 29,342 have collocated satellite data. 

The visual data were normalized and converted to albedos 
using the Muench and Keegan (1979) normalization scheme. 
The data set was separated into four categories 
(precipitation/no-precipitation, and inf rared/visuai) and 
after testing for normality, it was determined that none of 
the categories were normally distributed. Using histograms, 
a distinct separation between the peaks of precipitating and 
non- precipitating events was found, but some overlap does 
exist . 

Testing of infrared/visual thresholds for precipi taion/ 
no-precipitation events used in automated cloud and precipi- 
tation research yielded a correct estimation rate of 92% 
when the infrared and visual thresholds were combined. 
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ISIhodoc^on 

Every day, military aad civilian forecasters are tasked 
to provide weather forecasts for remote areas of the world. 
Many are short-range forecasts (up to 12 h) for specific 
locations, rather than large areas. Facsimile charts 
combined with climatology provide reasonable wind and temp- 
erature forecasts, which can be cross-checked using satel- 
lite imagery. In addition, satellite imagery aids in 
preparing sky coverage forecasts. A major remaining problem 
is precipitation. Since there may be few, if any, timely 
surface observations available, it can be difficult to 
determine whether precipitation is falling beneath areas of 
cloud cover- Hhat is required are methods to differentiate 
precipitation from non-precipitation using geostationary 
visual and infrared satellite data. 

Research in this area is important for two inter- 
changeable reasons. First, better observations will lead to 
better forecasts. A fundamental problem with all fore- 
casting schemes is a lack of reliable initial data. An 
accurate and reliable scheme for differentiating between 
precipitating and non- precipitating clouds will greatly 
improve current and future forecasting schemes. Second, the 
use of satellite data will increase the area of coverage and 
therefore increase the number of potential observations. 
The more "good” initial data the better the chance the fore- 
casting scheme will have to produce a better forecast. 

Although microwave sensing has shown much promise in 
detecting rainfall, visible and infrared data are utilized 
in this research. A major drawback of microwave sensing is 
that it uses a broad field of view and has a more coarse 
resolution. In addition, visible and infrared data are 
currently more readily available (every 30 minutes) . 
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Previously, studies have been undertaken which have 
attempted to correlate precipitation with satellite imagery. 
Muench and Keegan (1979) studied GOES-E visual and infrared 
satellite imagery and hourly rainfall climatological data 
for five surface stations in the northeastern United States 
for the period April through Novemoer 1977. Iheir data set 
consisted of 2760 surface observations, of which 300 
reported precipitation. Using infrared and visual imagery, 
they attempted to predict hourly rainfall amounts. Although 
their results were "far from ideal", they were still "useful 
for making decisions if the threshold for 'action' versus 
•no action' was known" (Huench and Keegan, 1979) . 

In addition to Huench and Keegan (1979), other 
researchers have studied this problem. Liljas (1981a, 
1981b) developed a bi-spectral cloud classification based on 
visual and infrared data from the polar-orbiting TIBOS-6 
satellite. Lovejoy and Austin (1979), using the results of 
Muench and Keegan (1979), developed a two dimensional 
frequency plot for cumulus and non-cumulus precipitation. 
In addition, they also tested a spectral threshold technique 
for rain area mapping. For those unfamiliar with research 
in this area, a summary by Paul (1983) of selected studies 
is included in Appendix A. 

Paul (1983), using only daylight observations (0800L - 
1600L EDI), computed means and standard deviations for each 
10 X 10, 8x8, 6 X 6 and 4 x 4 matrix of pixels (both 
infrared and visual) centered on a surface reporting station 
with collocated geostationary satellite data. Her results 
showed a distinct separation between the means of the 
precipitation and non-precipitation cases. However, adding 
one standard deviation to the infrared precipitation mean 
and subtracting one standard deviation from the non- 
precipitation mean resulted in a significant overlap between 
precipitation and non-precipitation thresholds. This 



overlap was too large to permit an unbiased forecast to be 
made. Similar results were found using visual satellite 
data . 

This research will concentrate on relating satellite 
imagery to observed precipitation and non-precipitation 
events using the surface-satellite data set of Paul(1983). 
The objectives are: 

1. Use a full 24 h data base, vice the 9 h data base by 
Paul (1 983) ; 

2. form distributions of satellite data for several 
classes of observed weather events, test for normality and 
study the separation of the distributions. 

3. Apply established visual and infrared thresholds to 
the data set and determine their reliability. 

This thesis is organized into four chapters. Chapter II 
describes the data set. Chapter III presents an analysis 
and discussion of the data set. Selection of matrix sizes 
for analyses and choices for studying the characteristics of 
each distribution are explained. Results are then compared 
to Paul (1983). The final section verifies the 
precipitation/no-precipitation infrared and visible thresh- 
olds currently used in SPADS. Chapter IV states conclusions 
and suggests further study. 
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II. DATA DESCEIPTION 



The data set assembled for this study consists of collo- 
cated GOES-E satellite data and Service-A hourly surface 
observations for the southeastern United States during 
August 1977. During August/ this region is dominated by 
subtropical airmasses with extensive convective activity. 
The GOSS-E data consists of 10 x 10 pixel matrices of 
visible and infrared satellite data centered over each of 
137 reporting stations (Fig, 1) , all south of 40 N. The 
satellite data are measured with the Visual Infrared Spin 
Scanned Eadiometer (VISSE) which have subsatellite point 
spatial resolutions of 1 km for visual channels and 7 km for 
infrared channels (U. S, Department of Commerce, 1983) The 
GOES-E navigation was completed by Man-computer Interactive 
Data Access System (McIDAS) at the University of IJisconsin 
using the full resolution visual data, with an accuracy of 
1-2 pixels (1-2 km) . The full resolution visual data were 
averaged to a 7 km resolution, to egual the infrared data 
resolution. The visual and infrared digital counts ranged 
from values of 0-63 and 0-255, respectively. The 10 x 10 
GOES-E visual and infrared satellite data each cover an area 
of 45 n mi X 45 n mi at 30 N (60 n mi x 60 n mi at 42 N). 

What distinguishes this data set from those used in most 
other studies is the number of observations. This data set 
has a total of 70,623 Service-A hourly reports, of which 
29,342 reports have collocated visual and/or infrared satel- 
lite imagery. The visual data were normalized and converted 
to albedos based on the algorithm of Muench and Keegan 
(1979). Their scheme corrects for the varying zenith angles 
as well as adjusting the visual satellite data for aniso- 
tropic scattering as related to the zenith angle. While 
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this scheme allows values 



albedo values cannot exceed 1.00, 
to overshoot 1.00 up to a value of 1.20. Therefore, the 
visual satellite values are not true albedos, but are esti- 
mated albedos. The extended visual normalized data scale 
was used to permit comparison of the results of this 
research to the bi-spectral threshold specifications of 
Muench and Keegan (1979) and to Paul (1983). The Muench and 
Keegan normalization scheme specifies that any computed 
albedo greater than 1.20 be set equal to 1.20 to limit 
unreasonably large values. Similarly, the scheme specifies 
computed albedos less than 0.15 be interpreted as the ground 
or water surface reflectance and the value 0.00 be assigned. 
(See Appendix B for further specific information concerning 
the Muench and Keegan normalization scheme.) The infrared 
data were processed in digital counts and converted to cloud 
top temperatures prior to statistical calculations. 

The surface reports within the data set contain informa- 
tion that regularly appears in ail Service-A reports, plus 
additional location information and visual and infrared 
satellite data (see Table I) . Of the information available, 
we were most concerned with the current weather, the cloud 
group, and the satellite data. 

For this study, the data set was divided into four 
categories: 

1. Observations with no precipitation and infrared 
satellite data (19,354 Obs) ; 

2. Observations with no precipitation and visual satel- 
lite data (14, 740 Obs) ; 

3. Observations with precipitation and infrared satel- 
lite data (1856 Obs) ; and 

4. Observations with precipitation and visual satellite 
data (1437 Obs) . 

The no-precipitation data category is comprised of all 
observations not showing any "E" in the current weather 
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group. The precipitation data category consists of all 
observations with an ”R” in the current weather group. 

The cloud cover group is based on a three digit code. 
The first digit indicates the amount of low clouds, where 0 
is defined as clear, 1 is scattered, 2 is broken, and 3 is 
overcast. The definitions for clear, scattered, broken, and 
overcast are as defined in the Federal Meteorological 
Handbook No. 1 (U. S. Department of Commerce, 1980) . The 
second and third digits of the cloud group indicate the 
amount of middle and high clouds respectively. 

In this study there were no restrictions to a specific 
time frame. Rather, all observation, regardless of time, 
which had satellite data available were used. The only 
restrictions imposed were: 

1. Only broken or overcast skies were considered; and 

2. Each pixel count must be greater than zero. 

Only broken or overcast ceilings were chosen to be analyzed, 
based on the assumption that precipitation from scattered 
skies is rare and that precipitation from clear skies is an 
incorrect observation. Indeed, Paul (1983) reported that 
less than one percent of the precipitation observations 
within her time period also reported scattered skies. The 
restriction that each pixel count be greater than zero was 
imposed because any pixel less than zero in meaningless. 

Paul (1983) limited her time period to daylight hours 
only (0800 - 1600L EDT) to avoid distortion of the visible 
data due to low solar elevation angles. By choosing not to 
limit the time period for this study, it was understood that 
distortion of visible data may occur. 
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III. DATA ANAkISIS 



A. INTRODOCTION 

This chapter consists of analyses and discussions of the 
data set. It begins with a discussion on choosing the most 
representative and reliable matrix size and the distribu- 
tions formed by those matrices. Next, the distributions are 
tested for normality. Using those results, the distribution 
characteristics are described and discussed. Next the 
results are compared to those obtained by Paul (1983). In 
the final section, the data set is used to verify the estab- 
lished precipitation/no-precipitation thresholds of an auto- 
mated cloud and precipitation estimation scheme (Hash eX al , 
1984) with the U. 5. Navy’s Satellite Data Processing and 
Display System (SPADS). 

B. MATRIX SIZE 

The collocated satellite data consisted of a 10 x 10 
matrix centered on the surface observation station. This 
provided the option of using any one of nine different 
matrix sizes, 10 x 10 through 2x2, in the research. Paul 
(1983) evaluated four sizes, 10 x 10, 8x8, 6x6 and 4x4 
(2025 sg n mi, 1296 sg n mi, 729 sg n mi and 484 sq n mi 
respectively). Her results (using 7358 no-precipitation 
observations and 534 precipitation observations) showed no 
significant differences in the means and standard deviations 
of the four matrices. Hhen Paul's time restrictions were 
removed, the results (using up to 19,354 no-precipitation 
observations and 1856 precipitation observations) also 
showed no significant differences between the same four 
matrices based on the means and standard deviations (Table 
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II). From these results it was concluded that the distribu- 
tion of the data would not change significantly as the 
matrix sizes were changed. The matrix size could now be 
based on meteorological parameters. 

Initially one could argue that the 2x2 matrix would be 
the most representative. It could most accurately represent 
the surface observation by not being influenced as much by 
the surrounding area as any of the larger sizes. The major 
disadvantage of the 2x2 matrix though, is actually its 
small size. Because the navigational error of dcIDAS is one 
to two pixels, the possibility that the surface observation 
would actually be outside the satellite data matrix is 
present. Choosing the 4x4 matrix would greatly reduce the 
possibility of having the surface observation outside the 
satellite data matrix. The 6 x 6 matrix would guarantee 
(within McIDAS error) the surface observation being within 
the matrix, but now areal coverage must be considered. The 
6x6 matrix has an areal coverage of 729 sg n mi at 30 N 
while the 4x4 matrix covers 484 sg n mi, approximately 66% 
of the 6 X 6 matrix. After weighing the possibility of a 
misplaced observation using the 4 x 4 matrix, with the 
increased coverage of the 6x6 matrix, the 4x4 matrix was 
chosen as the most representative matrix size. 



C. NORMALITY 

The determination of the normality of the distributions 
within the data set was an important stipulation of this 
research. Since the normality determined the way in which 
the distributions could be described (using confidence 
levels if normal, and just general descriptions if not 
normal), it was paramount that the distribution be tested 
correctly. 
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Histograms were plotted for eacL of the four data 
categories with observations wnich had collocated satellite 
data. Within each category, a 10 x 10, 8 x 8, 6 x 6 and 4 x 
4 matrix was also plotted. Even though the 4x4 matrix had 
already been chosen as the size to be evaluated, the option 
to change sizes if warranted by the normality testing was 
left open. 

The Chi-sguare test was chosen to test the distribution 
for normality. It involves comparing a sample distribution 
to a normal distribution. Each distribution was divided 
into 20 bins. Using the mean and standard deviation 
computed earlier for each category, 17 degrees of freedom 
were available. All of the no-precipitation categories 
could be seen not to be normally distributed by inspection 
of the histograms (see Fig. 2 and Fig. 3) and therefore were 
not formally tested. Since each precipitation category more 
closely resembled a normal distribution, they were tested 
with the Chi-sguare test. Each no-precipitation distribu- 
tion tested (all four matrix sizes) proved not to be 
distributed normally. Fig. 4, although appearing to be 
normally distributed, only had a 5% chance of coming from a 
normal distribution, and it was the closest to normal of any 
category tested. 

It must be understood that just because a distribution 
has a low confidence level for being normally distributed, 
it does not mean that it is not a normal distribution. 
Likewise, just because a distribution has a high confidence 
level of being normally distributed, does not mean that it 
is a normal distribution. Since the results showed a high 
probability of not being normally distributed, the research 
focus shifted to describing the non-normal and other inter- 
esting aspects of the distributions. 
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D. DISTEIBOTIOH CHABACTEBISTICS 



1 . Introdu cti on 

This section will present a dxscussion about each of 
the histograms plotted for xnfrared/v isu al and 
no-precipitation/precipitation. It is important to remember 
that only broken or overcast ceilings are included in the 
distributions. Also, due to round-off error, not all 
percentages calculated in this section add up to 100% 
exactly. 

2- No-Pr ecipi ta tio n, Infrared 

This category consists of all surface observations 
without precipitation and with collocated infrared satellite 
data (see Fig. 2). There were 19,354 observations within 
this category. As expected, the majority of the observa- 
tions are grouped near the warm end. .The cloud top tempera- 
tures ranged from a minimum of 205.0 K to a maximum of 320.1 
K. It has a peak of 4214 observations between 290.0 K and 
295.0 K. The mean temperature was 275.7 K with a standard 
deviation of 19.4 K. 

The most interesting part of this distribution is 
the cold tail (less than 240. 0 K) . These are observations 
which would be misclassif ied by a single cold threshold. It 
was found that of the 1461 observations within the cold 
tail, overcast and broken ceilings were divided evenly (see 
Table III). Within the overcast ceiling category, 35% were 
high clouds with 10% mid- level clouds and 5% low clouds. 
Within the broken ceiling category, 41% were high clouds 
with 7% mid-level clouds and 3% low clouds. If two or more 
levels reported broken ceilings, the highest ceiling was 
specified as the ceiling for this study. For example, if a 
station reported 3500 broken, 10,000 broken, and 22, 000 
broken, the ceiling would be reported as a high ceiling. 
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Combining the overcast ceilings with the broken ceilings 
produced 76% of the observations with high clouds, 16% with 
mid-level clouds, and 8% with low clouds. 

The breakdown of the overcast ceiling category can 
be misleading. When an observer reports an overcast 
ceiling, he cannot see what is above the base of the clouds. 
The misrepresentation arises because the tops of the clouds 
are not at the level reported by the observer. In addition, 
there may be additional layers above the overcast layer. 
The satellite data will therefore contain information about 
the highest (coldest) clouds it senses and that layer may 
possibly be one or two layers above the reported ceiling. 

To approximate the percentage of overcast ceilings 
affected, the high cloud with broken ceilings category was 
further analyzed. It was found that 95% of the observations 
within that category were multi-layered clouds (with either 
scattered or broken conditions reported at low or mid-levels 
in addition to a broken ceiling reported at the high level) . 

The overwhelming number of warm cloud top tempera- 
tures is not surprising for a non-precipitation category. 
Of some concern though, are the few temperatures (66) which 
exceed 300.0 K (27.0 C). One explanation is that the 
infrared sensor is seeing the top of a marine layer along 
the Gulf coast. Another possibility is that the surface of 
the earth is influencing the mean value of the 4x4 matrix. 
This is possible because one surface observation is repre- 
senting a 424 sg n mi area. Just because a ceiling is being 
reported at the observation location, does not necessarily 
mean there is a ceiling 10 n mi away. Warm surface tempera- 
tures could increase the mean value of the matrix 
significantly. 

In summary, the large majority of warm temperatures 
for non-precipitation cloud tops is just as one would 
expect. The extremely warm temperatures do not affect the 
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precipitation/no-precipitation decision. The cold tail is 
composed mostly of thick hign clouds or overcast/broken 
clouds with additional clouds at lower levels. Significant 
errors arise if these cold cloud tops are interpreted as 
precipitating clouds. 

No-Pr ec ipita tion^ Visual 

This category consists of all observations without 
precipitation and with collocated visual satellite data (see 
Fig. 3). There were 14,740 observations within this 
category. As expected, the majority of the observations are 
grouped near the dim end. The estimated albedos ranged from 
0.00 to 1.20. It has a peak of 1944 observations at 0.00. 
The mean albedo was 0.27 with a standard deviation of 0.26. 

Further analyses were conducted on the 999 observa- 
tions comprising the bright tail (greater than 0.70) which 
would lead to misclassif ication of precipitation from the 
visual data. It was found that 49% of the observations 
reported overcast ceilings and 51% reported broken ceilings 
(see Table IV) . Within the overcast ceiling category, 18% 
were high clouds, 14% were mid-level clouds, and 17% were 
low clouds. Within the broken category, 41% were high 
clouds, 7% were mid-level clouds, and 3% were low clouds. 
Combining the overcast ceilings with the broken ceilings 
produced 59% of the observations with high clouds, 22% with 
mid-level clouds, and 20% with low clouds. As in section 
D.2, the breakdown of the overcast ceilings can be 
misleading. To approximate the actual percentage of high 
clouds, further analyses were conducted on the high cloud 
with broken ceiling category for multi-layered clouds. It 
was found that 83% of those observations contained multi- 
layered clouds. 

The large number of albedos between 1.15 and 1.20 is 
cause for concern. Paul (1983), by imposing a time 
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restriction on the data base. excluded nost of the problems 
associated with distortion caused by low solar angles near 
sunrise and sunset. Because this study did not restrict the 
data base to a specific time period, distortion was a 
problem. In the distribution for this category there were 
336 reports with albedos greater than 1.15. Further anal- 
ysis of these reports showed all albedos greater than 1.15 
were reported at 0000 GMT (near sunset) . In addition, 83% 
of all observations with albedos greater than 0.70 occurred 
within three hours of sunrise or sunset. In addition to 
high albedo values caused by distortion or nearby precipita- 
tion, further analysis showed that one- third (44/135) of the 
observations betweem 0.70 and 1.15 which reported low over- 
cast ceilings also reported fog. It must be noted that 
although these observations were clusrered around sunrise, 
they were still assigned higher albedo values than one would 
normally expect for fog. A significant portion of these 
bright non-precipitation observations are related to prob- 
lems in normalizing the satellite data for low sun angle 
situations. More work is reguired in using the cloud 
brightness normalization before more quantitative use of the 
data can be made. 

In summary, most of the bright tail is composed of 
thick high clouds or overcast/broken clouds with clouds also 
at lower levels. As in the no-precipitation, infrared 
category, significant errors arise when these bright clouds 
are interpreted as precipitating clouds. In addition, 
distortion due to low solar angles also causes misclassif ied 
bright albedo values . 

Precip it atio n, Infrare d 

This category consists of all observations reporting 
precipitation and with collocated infrared satellite data 
(see Fig. 5) . • There were 1856 observations within this 
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category. The distribution of cloud top temperatures (clus- 
tered near the cold end) is similar to what would be 
expected for a precipitation category. The cloud top temp- 
eratures ranged from a minimum of 205.0 K to a maximum of 
296.6 K. It has a peak of 181 observations between 220.0 K 
and 225.0 K. The mean is 243.0 K with a standard deviation 
of 21.7 K. 

Further analyses were conducted on the warm tail 
(warmer than 270.0 K) , the cold end (colder than 240.0 K) , 
and the area between them (see Table V). The warm tail was 
analyzed because it would most likely contain misclassif ied 
observations. The cold end and the middle area each were 
analyzed to better understand the satellite data distribu- 
tion. Each of the three areas were analyzed for the type 
and intensity of the precipitation, and for convective 
versus stratiform precipitation. 

In the warm tail there were 257 observations. Light 
rain showers were dominate with 50% of the observations. 
Light rain and light thunderstorms followed in frequency 
with 21% and 18%, respectively. An overwhelming majority, 
91%, reported light precipitation, compared to only 9% with 
moderate/heavy precipitation. In addition, 76% of the 
observations were convective in nature compared to 24% of a 
stratiform nature. One can expect the observations within 
the warm tail to be misclassif ied because precipitation is 
not usually associated with warm cloud top temperatures. It 
should be noted though, that only 9% of the precipitation 
observations within the warm tail reported moderate or heavy 
precipitation, whereas 91% reported light precipitation. 
This leads one to conclude that although the observations 
within the warm tail could be misclassified, they would 
would be , light precipitation and not moderate/heavy 
obs ervat ions. 
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some 



A possible reason why some precipitation 
precipitation observations have such warm cloud tops is 
because the precipitation, in addition to being light, is 
also very scattered. In this situation, a portion of the 
satellite data is composed of surface reports rather than 
cloud top data. Also, it should be noted that the surface 
observations are taken 5-10 minutes before the hour, while 
the satellite scan is taken 5 minutes after the hour. This 
10 to 15 minute lag in the satellite scan is at least 
partially responsible for some of the misclassif ied observa- 
tions. Precipitation could be reported at the surface 
station at the observation time, but by the satellite scan 
time, it may have stopped. This would produce a surface 
precipitation report with a collocated no-precipitation 
satellite report. 

In the cold end there were 959 observations. light 
thunderstorms were reported 42% of the time with light rain 
showers being reported 20% of the time. In addition, 15% of 
the observations reported moderate/heavy thunderstorms and 
14% reported light rain. There were more light precipita- 
tion reports (81%) than moderate/heavy reports (19%). Once 
again convective type precipitation (79%) was reported more 
often than stratiform precipitation (21%) . 

In the region between the warm tail and the cold end 
there were 640 observations. Light rain showers were in 37% 
of the observations, followed by light rain in 29% and light 
thunderstorms in 21% of the observations. Light precipita- 
tion was reported in 89% of the observations, with moderate/ 
heavy precipitation in the remaining 11%. As before, 
convective type precipitation (66%) was reported more often 
than stratiform precipitation (34%) . 

It should be noted here that this study cannot 
distinguish between precipitation intensities with any 
degree of accuracy. One would expect a much higher 
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percentage of moderate/heavy thunderstoras in the cold end 
than in the wara tail. These results showed only a slight 
increase [ 9 % to 19%) in occurrence of heavier precipitation 
moving from the warm tail to the cold end. The major reason 
is that this study uses an average of surrounding points to 
generate a collocated satellite observation, whereas other 
studies, which show a distinct increase in precipitation 
intensity with colder cloud top temperatures, used the 
coldest cloud top temperature as the value. 

In summary, although the occurrences of moderate/ 
heavy precipitation were less than expected, the trend of 
increasing intensity with colder cloud top temperatures was 
present. Also, the trend of increasing thunderstorms and 
decreasing rain showers with colder cloud top temperatures 
is consistent with expectations. 

5. Preci pit at ion, Visi ble 

This category consists of all observations which 
reported precipitation and had collocated visible satellite 
data (see Fig. 6). There were 1437 observations within this 
category. The distribution is as one would expect for a 
precipitation with visible satellite data category (albedo 
values clustered near the bright end) . The estimated 
albedos ranged from 0.00 to 1.20. There is a peak of 132 
observations between 0.65 and 0.70. The mean albedo is 0.61 
with a standard deviation of 0.29. 

This distribution was divided into three areas and 
further analyses were conducted on the dim tail (albedos 
less than 0.40) , the bright tail (albedos greater than 
0.90), and the region between the dim and the bright tails 
(see Table VI) . The dim tail was chosen because it is the 
area of the greatest probability of misclassif ied observa- 
tions. The other two areas were chosen to compare to the 
din tail. Within each of the three areas, an analysis for 
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precipitation type, precipitation intensity and convection 
or stratiform precipitation was conducted. 

In the dim tail there were 327 observations. Light 
rain showers (43%) were dominate with light thunderstorms 
(28%) and light rain (15%) following in frequency. Light 
precipitation was reported in 88% of the observations with 
moderate/heavy precipitation in 12% of the observations. In 
addition, convective type precipitation was reported 81% of 
the time with stratiform precipitation being reported 19% of 
the time. As expected, most of the misclassif ied observa- 
tions were light precipitation reports. 

In the bright tail there were 165 observations. 
Light thunderstorms were reported in 29% of the observations 
with light rain being reported in 21% of the observations. 
Also, 19% of the observations were light rain showers and 
18% were moderate/heavy thunderstorms. Light precipitation 
was reported more often (72%) than moderate/heavy precipi- 
tation (28%) . Convective type precipitation was also 
reported in 72% of the observations as opposed to stratiform 
precipitation being reported in 28% of the observations. As 
expected, most of the moderate/heavy precipitation and thun- 
derstorm observations were in the bright end. 

In the region between the dim tail and the bright 
tail, there were 947 observations. Light rain showers and 
light thunderstorms were dominate in this region with 30% 
and 29% of the reports respectively. Light rain was 
reported in 22% of the observations. Light precipitation 
was reported in 83% of the observations while moderate/heavy 
precipitation was reported in 17% of the observations. 
Within this category, convective type precipitation was 
reported 74% of the time while stratiform precipitation was 
reported 26% of the time. 

An important consideration when evaluating the 
visual categories is the time of day. In this category, 57% 
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of the observations were within a tnree hour period from 
2200 to 0000 GMT. This implies that most of the 
precipitation occurs late in the afternoon. Also, 37% of 
the 0.00 estimated albedo values (39) occurred at 0000 GMT. 
Of the remaining 0.00 values, two occurred at 1300 GMT and 
one each at 2000 GMT, 2100 GMT, and 2300 GMT. This suggests 
that distortion caused some of the very low estimated albedo 
values. 

E. CCMPABISON OF EESOLTS TO PAUL (1983) 

1 . Intro du cti on 

This section will compare our results with the 
results by Paul (1983). Paul (1983) divided her data into 
seven precipitation categories: 

1. convective; 

2. continuous; 

3. light; 

4. moderate/heavy; 

5. general; 

6. overcast ceiling; and 

7. overcast and broken ceiling; 
and two no-precipitation categories: 

1. no-precipitation, overcast ceiling; and 

2. no-precipitation, broken and overcast ceiling. 
Each of the above nine categories were further divided into 
those with infrared satellite data and visible satellite 
data . 

As stated previously, this study used one precipita- 
tion category and one no-precipitation category, each 
divided into a visible and an infrared satellite group. In 
addition, only those surface observations which reported 
either broken or overcast ceilings were evaluated. 
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This chapter is divided into two sections. One will 
compare results of the precipitation category and the other 
will compare results of the no-precipitation category. 

2. Pr e cip i tat ion 

This section will compare these results to those 
from Paul (1983) for the precipitation with broken or over- 
cast ceilings (see Table VII) . The mean values will be 
discussed first, followed by the standard deviations. 

a. Means 

The mean infrared values of this study are 
approximately 8.0 C colder than those from Paul (1983) (see 
Table VII) . The differences can be attributed to the fact 
that nightiae data are included in this study. In this area 
of the country many thunderstorms occur from late evening 
through the early morning hours. By including these storms 
(with their cold cloud top temperatures) in this study, the 
mean cloud top temperatures must decrease. 

Also, the means of the infrared data decreased 
slightly (became colder) in both studies as the matrix size 
decreased. This decrease can be attributed to the decrease 
in surrounding area when changing from a 10 x 10 matrix to a 
4x4 matrix. The smaller matrix size allows the satellite 
to detect only the colder cloud tops associated with the 
precipitation and not the warmer clouds from the fringe 
areas . 

There was very little difference in the visible 
mean values between the two studies (see Table VII) . This 
is not surprising because the time restrictions imposed by 
Paul (1983) contained most of the visible satellite data 
within the data set. The trend of the mean visual values is 
the same for both studies (see Table VII) . Both show a 
slight increase in albedo values when changing from a 10 x 
10 matrix to a 4x4 matrix. 
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b. Standard Deviations 



The infrared standard deviations are slightly 
larger in this study tnan in Paul (1933) (see Table VII) . 
The trend of increasing standard deviations with decreasing 
matrix size occurs in both studies. This can be related to 
the warmer fringe areas of the clouds being deleted by the 
smaller matrices. 

The standard deviations of the albedo values are 
larger in this study than in Paul (1983) (see Table VII). 
By restricting the time frame, Paul (1983) eliminated many 
possible distorted values that were included in the calcula- 
tions for this study. Since the distorted values tend to 
accumulate at the extremes (0.00 and 1.20), including them 
in the calculations increases the standard deviations. 

The trend of the albedo standard deviation 
values is toward slightly larger values as the matrix sizes 
decrease. The less fringe area covered by the 4x4 matrix 
diminishes the chances of the dimmer fringe area values 
being included in the calculations. 

3- No-precipita tio n 
a. Means 



The mean infrared values in this study are 
approximately 4,0 C colder than those from Paul (1983) (see 
Table VIII) . The mean values of both studies remain almost 
constant as the matrix size is decreased. 

The visual means showed very little difference 
between the two studies (see Table VIII). As stated before, 
this is not surprising since most of the visual readings are 
within the time restrictions imposed by Paul (1983). In 
both studies, the mean albedo values tend to remain constant 
as the matrix size is decreased. 
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b. Standard Deviations 



The infrared standard deviations are about 1.5 C 
larger than Paul ( 1983) (see Table VIII) . As in the 
precipitation category the standard deviations increase 
slightly as the matrix sizes decrease. 

The standard deviations of the visible values 
are approximately 0.06 larger than those in Paul (1983) (see 
Table VIII) . This can be attributed to the addition of 
distorted mean values in this study which were excluded by 
the time restrictions imposed by Paul (1983). The distorted 
values give additional weighting to each end of the range, 
causing the standard deviations to increase. The trend of 
the slightly increasing visual standard deviations with 
decreasing matrix sizes is again repeated here. 

F. SPADS THEESHOLD VERIFICATION 

This section will verify the no-precipitation/ 
precipitation threshold of a cloud and precipitation anal- 
ysis model of the United States Navy's interactive Satellite 
Data Processing and Display System (SPADS) using the 24 h 
data base. The objective is to present a comprehensive 
analysis of the infrared, visible, and infrared and visible 
thresholds combined. No attempt will be made to derive 
better thresholds, but only to verify the existing thresh- 
olds. After a brief explanation of the automated cloud and 
precipitation analysis program, there is a discussion on how 
thresholds work in an ideal situation. Next is an explana- 
tion on how the verification was conducted, and then the 
results are presented. 

This program was designed to produce in real time (15 to 
30 minutes) analyses of important cloud and weather 
features. It uses the Geostationary Observational 
Environmental Satellite (GOES) visual and infrared images to 
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produce contoured digital displays of cloud types, cloud 
amounts, cloud top heights, cloud top temperatures, and 
precipitation intensity for an approximate 1024 x 1024 n mi 
sector of a geostationary image. 

The objective was to determine the success of the 
precipitation/no-precipitation thresholds in separating 
precipitation from no-precipitation for this data set. 
Ideally a threshold value will produce a clear and distinct 
separation between two events. If the two distributions 
overlap (an area with egual probability of either event 
occuring) , the smaller the overlap the better. The 
threshold with an overlap can still be useful as long as 
there is not too high a percentage of each event located 
within the overlap region. The more events located within 
the overlap region increases the chance of a misclassif ied 
event up to a point where there is no longer any skill asso- 
ciated with this prediction technique. 

The verification was conducted using the entire data set 
over the full 24 hours. The visual and infrared thresholds 
from the interactive cloud and precipitation analysis 
program were verified using the four categories described in 
chapter II. The thresholds values used in the verification 
were : 

1. Estimated albedos greater than or equal to 0.55; and 

2. Cloud top temperatures less than or equal to 249.0 K; 
for precipitation, and: 

1. Estimated albedos less than 0.55; and 

2. Cloud top temperatures greater than 249.0 K; 
for no-precipitation. 

Using only the infrared threshold, the no-precipitation 
and precipitation histograms (see Fig. 2 and Fig. 5) were 
scanned for values which exceeded/were less than or equal to 
249.0 K. The results (see Table IX) showed that 64% of the 
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precipitation observations had cloud top temperatures 
greater than or equal to 249.0 K. Likewise, 87% of the 
no-precipitation observations had cloud top temperatures 
less than 249.0 K. 

Using only the visual thresholds, the no-precipitation 
and precipitation histograms (see Fig. 3 and Fig. 6) were 
scanned for values which were greater than or equal to/less 
than 0.55. The results (see Table X) showed that 63% of the 
precipitation observations had albedo values greater than 
0.55. Likewise, 87% of the no-precipitation observations 
had albedo values less than 0.55. 

Next the above calculations were repeated, only this 
time the criterion for precipitation/no-precipitation was 
based on both the visual and infrared satellite data. If an 
observation did not contain both visual and infrared satel- 
lite data, it was not considered in our calculations. The 
results are shown in Table XI. They show that 82% of all 
precipitation observations in the data set had both albedo 
values greater than or equal to 0.55 and cloud top tempera- 
tures less than or equal to 249.0 K. Likewise, 93% of the 
no-precipitation observations had albedos less than 0.55 and 
cloud top temperatures greater than 249. 0 K. Using the 
visual/infrared thresholds together increases the chance of 
correctly classifying a precipitation event from 63% for 
visible and 64% for infrared to 82%. In addition, using the 
visual/infrared thresholds together increases the chance of 
correctly classifying a non-precipitation event from 87% for 
both visual and infrared to 93%. 

G. SUMMARY 

This chapter consisted of analyses and discussions of 
our data set. It began with a discussion on choosing the 
most representative and reliable matrix size and the 
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Next the distribu- 



distr ibutions formed by those matrices, 
tions were tested for normality. Using those results, a 
discussion of the characteristics of the distributions were 
presented. From there, the results were compared to those 
obtained by Paul (1983) . The final section was a verifica- 
tion of the established precipitation/no-precipitation 
thresholds using the data set. 
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IV. SDMHAfiy AND CONCLOSIONS 



A. DATA DESCRIPTION SOMHAEY 

A data set consisting of 70,623 surface observations, of 
which 29,342 had collocated satellite data, was studied. 
The satellite data consists of a 10 x 10 matrix centered on 
the surface observation station. The visual satellite data 
were converted to albedo values using the Huench and Keegan 
(1979) normalization scheme. This scheme also corrects for 
anisotropic scattering effects. The data set was divided 
into four categories of precipitation/no-precipitation and 
visual/infrared satellite data. It was further divided by 
using only those observations which reported broken or over- 
cast ceilings. In addition, the full 24 hours of available 
data was analyzed, compared to the 9 hours of daytime obser- 
vations used by Paul (1983) . 

B. DATA ANALYSIS SOMHARY 

This chapter consisted of analyses and discussions of 
the data set. It was determined that the 4x4 matrix had 
the best combination of accuracy of location, and accuracy 
in best describing the region based on a single observation. 
Next, a description of how the histograms were formed and 
tested for normality was given. Because one could not say 
with any acceptable degree of confidence that the distribu- 
tions were normally distributed, it was decided to analyze 
the non-normal parts and to gain a better understanding into 
the composition of the data set. The trends when moving 
from the cold/bright end to the warm/dim ends were as 
expected, but one could not accurately separate light from 
moderate/heavy precipitation using the cloud top 
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temperatures or albedo&2)When the results were compared to 
Paul (1983), there were many similarities, especially when 
decreasing the size of the matrices from 10 x 10 down to 4 x 

4. 

The last area studied was the verification of estab- 
lished precipitation/no-precipitation thresholds used by the 
cloud and precipitation analysis scheme in SPADS with the 
data set. It was found that using both the visual and 
infrared satellite data together to determine precipitation/ 
no-precipitation produced the best results. 

C. CONCIOSIONS 

The conclusions of this study are: 

1. Based on the results of the Chi-sguare test, the 
data set and all its subsets used in this study are not 
normally distributed; 

2. The histograms for the precipitation/no- 
prec ipit ation and infrared/visual categories each were 
distributed as one would expect for precipitation and 
no-precipitation events respectively; 

3. Although the peaks of the precipitation/no- 
precipitation and infrared/visual distributions were as 
expected, the precipitation and no-precipitation distribu- 
tions display partial overlap; 

4. Investigation of the overlap areas indicates that 
multiple and thick high cloud layers are responsible for 
erroneous estimates of precipitation within the 
no-precipitation distributions. Concerning the precipita- 
tion distributions, the scattered character of the summer- 
time precipitation likely produced a significant number of 
misclassifications ; 

5. Albedo estimates for the brightness normalization 
scheme are not reliable in the low sun angle situations; and 
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6. Estimates of reliability for precipitation/no- 
precipitation thresholds were obtained from the data set. 
Using both the infrared and visual precipitation/no- 
precipitation thresholds combined greatly improves the 
precipitation estimation when using each threshold sepa- 
rately. 

D. SUGGESTED FUBTHEE STUDY 

The recommendations for further study are: 

1. Further analyses of precipitation categories should 
be conducted by separating them into specific precipitation 
types (i-e- thunderstorms, rain showers, etc) ; 

2. Do additional studies of precipitation/no- 
precipitation thresholds for cloud/precipitation estimation; 
and 

3. Gather and study additional collocated data sets in 
other seasons and for other regions. 
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4££MMX k 

PREVIOUS RESEARCH BY PAUL (1983) 

A. BI-SPECTRAL AND INFRARED THRESHOLD 

Liljas (1981a, 1981b) developed a bi-spectral cloud 

classif ication based on visual and infrared data from the 
polar orbiting TIROS-6 satellite (see Fig. 7) . The data set 
consisted of a limited number of daily observations, chosen 
for their synoptic characteristics, in May and August 1979 
over a region encompassing Norway, Sweden, Finland, and the 
Baltic Sea with weather charts providing the ground truth. 
Based upon the precipitation threshold results of Muench and 
Keegan (1979), liljas chose a cloud top temperature 
threshold of -12 C to -15 C to classify cumulonimbus and 
nimbostratus clouds. Starting with this cloud classifica- 
tion and the assumption that the highest and the densest 
clouds produce the maximum precipitation amount, Liljas 
suggested a qualitative precipitation scale based on the sum 
of the visual and infrared satellite digital counts (see 
Table XII) . These sums represent the areas of the Liljas 
nimbostratus and cumulonimbus cloud types in his bi-spectral 
cloud classification (see Fig. 8) . 

Lovejoy and Austin (1979) studied rain mapping of cloud 
areas based on GOES visual and infrared satellite data over 
Montreal, Canada, and the tropical Atlantic (Global 
Atmospheric Research Program Atlantic Tropical Experiment, 
GATE, data) with radar providing the ground truth. The 
Montreal data set consisted of 17 observations over three 
days during June 1977. Working with 4x4 resolution satel- 
lite image, Lovejoy and Austin plotted two dimensional 
frequency grids for the radar-determined rain and no-rain 
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points on a 25 X 25 array (See Fig. 9 and Fig. 10). The 
visual data were normalized by selecting the "brightest" and 
"dimmest" values in each image and linearly interpolating 
the radiances between 0 and 1. 

Lovegoy and Austin (1979) state, with reference to the 
cumulus rain data distribution of Fig. 9 that, "The distri- 
bution was to a good approximation a two-dimensional 
Guassian." They do not provide or describe the statistics 
to support this assertion. The no-rain cumulus cases (Fig. 
10) were described as a bimodal distribution with one peak 
near the lour visual and low infrared values and the other 
peak near the rain peak but shifted slightly toward lower 
values. In most cases, the separation of the cumulus rain 
and no-rain cases was statistically significant with the 
probability ranging from 10% to 50% that the rain and 
no-rain samples came from the same population. 

The Love joy and Austin (1979) two dimensional frequency 
plots for non-cumulus storms were limited to one case. The 
significant differences between the cumulus and non-cumulus 
data sets were the non-cumulus no-rain plot lost its bimodal 
character, relative to the cumulus no-rain plot, and 
appeared as a broad two dimensional Gaussian distribution. 
The non-cumulus rain plot points fell within the no-rain 
distribution, but were shifted slightly higher in the 
visual. The separation of the non-cumulus rain and no-rain 
cases was not statistically significant, with greater than a 
50% probability of the rain and no-rain samples coming from 
the same population. 

Love joy and Austin (1979) attempted to further classify 
the cumulus rain and no-rain cases into no-rain, light rain, 
and heavy rain. Rainfall rates greater than 2 mm-h-1, as 
determined by radar, were defined as heavy rain. As 
expected, the mean of the heavy rain cases was shifted 
slightly toward higher visual and infrared values than the 
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the shift was so 



dean of the light rain cases. However, 
saall that there was at least an 80% probability of the 
light rain and the heavy rain cases coming from the same 
population. Lovejoy and Austin (1979) concluded that 
•'little, if any rainfall-rate information is contained in a 
single (visual and infrared) satellite image. 

Lovejoy and Austin (1979) tested a spectral threshold 
technique for rain area mapping. Each satellite image of 
400 X 400 km was divided into one hundred 40 x 40 km boxes. 
The 100 subareas were each checked with radar to determine 
the total number of rain areas. An equal total number of 
satellite subareas were classified as rain areas. The 
satellite subareas with the highest visual and highest 
(cold) infrared values were classified as rain areas, until 
the total number of satellite rain areas equaled the total 
number of radar determined rain areas. This spectral 
threshold technique was applied to three days accumulation 
of data and is shown in Tables XIII and XIV. When compared 
to the success of the two dimensional frequency plot method, 
the visible and infrared thresholds averaged 45% and 58% 
worse, respectively. The accuracy of the visual threshold 
is limited by the extent of low, thick clouds and the 
infrared threshold is limited by the extent of the cirrus 
clouds in the satellite image. Lovejoy and Austin (1979) 
concluded that "the errors involved in using a ’best thresh- 
old' are very large indeed." 

Del Beato (1981) studied correlations between cloud top 
temperatures (based on NOAA-5 satellite data) and rainfall 
totals for 30 and 60 minute intervals over eastern 
Australia. The satellite data had a 60 sq km maximum reso- 
lution at subsatellite point and cloud top temperatures were 
area averaged for resolution of 200 sq km. The 21 data sets 
were first classified according to synoptic situation in a 
rough attempt to group the data by cloud type, droplet 
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spectra, and air mass trajectory. The initial results 
suggested that the cloud top temperature determined an upper 
limit on rainfall amount, with the maximum increasing as the 
cloud top temperature decreases. k linear correlation anal- 
ysis to determine a quantitative relationship between rain- 
fall amount and cloud top temperature gave indefinite 
results . 

Further study of surface and radiosonde observations 
indicated that classification by proportion of cumuliform 
cloud reports to all cloud reports and subcloud layer 
humidity might be more appropriate (Del Beato, 1981) . This 
classif ica tion resulted in a correlation coefficient of 
0.90, excluding cases with cumuliform portions less than 50T. 
and dew point depressions of greater than 6 C. Finally a 
composite frequency distribution was calculated based on 
three cases, all southwesterly stream situations described 
as ''post-frontal cellular convection cases in cyclonically 
curved flow." The fitted equation was: 

f = 0. 057 - 0. 004CT1 - 0.054E (A.1) 

where f is the rainfall frequency, R is the 30 minute rain 
total (tarn), and CTT is the cloud top temperature in degrees 
Celsius. The equation was fitted to 41 independent f 
values. This equation is associated with a correlation of 
0.79 at the 99% confidence level. Equation A.1 indicates no 
rain from clouds warmer than +13 C and a maximum 30 minute 
rainfall of 2.5 mm for a cloud top temperature of -20 C. 

In summary, Del Beato (1981) found that cloud top temp- 
eratures and 30 and 60 minute rainfall totals indicated 
statist ically significant relationships for cloud systems 
with a high portion of cumulus clouds and high subcloud 
humidity. Additionally, as cloud top temperatures decrease 
to at least -35 C, rainfall totals increase. 
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Wylie (1982) attempted to correlate rainfall occurrence 
with radiosonde soundings, hourly Service-A observations, 
and visual and infrared satellite data. His data sample was 
restricted to "large-scale cloud cover" areas with wide- 
spread precipitation (rain gauge reports varied less than 
20%) for the Great Plains States region for the period 27 
February 1981 through 4 January 1982. From 13 parameters 
derived from three data sources (see Table XV) the best 
linear regression equation for estimating rainfall rate was: 

6h rain = 1. 0242 + 0.380Pw - 0. 0304Qc - 0.0047Ct (A. 2) 

where pw is the vertically integrated precipitable water 
vapor (in) , Qc is the moisture convergence (g/kg/day) , and 
Ct is the cloud top temperature (Kelvin). Equation A. 2 has 
a linear correlation coefficient of 0.60. Linear regression 
equations were also determined for the three parameters 
alone and for a combination of Pw and Qc to be used when not 
all three data types were available. The cloud temperature 
regression equation was: 

6h rain (in) = 2. 10 - 0.008Ct (A. 3) 

The correlation coefficient was -0.35. V7ylie ( 1982) stated 
that the synoptic scale data base measurements were best 
suited for estimating broad changes in rainfall rates asso- 
ciated with air masses and not suited for estimating rain- 
fall rates associated with small scale dynamic processes. 

B. LIFE HISTORY 

The life history methods are empirically derived precip- 
itation estimation schemes based upon two assumptions, 
first, that significant rainfall comes from convective 



clouds, and second, that convective clouds can be identified 
and measured in satellite images. These methods involve 
manual analyses of convective cloud areas in a sequence of 
visual, infrared, or both visual and infrared satellite 
images. Threshold values and study condition parameters 
associated with published life history studies are summa- 
rized in Table XVI. 

The Scofield/Oliver (Scofield, 1981) analysis follows a 
decision tree procedure to estimate half-hourly rainfall for 
deep convective systems within tropical air masses. Using 
enhanced infrared and high resolution visual satellite 
images, the technique involves first identifying the active 
convective portion of the cloud, or cluster, from two 
consecutive satellite images. Once the active portion is 
identified, the half-hourly rainfall estimation is computed 
based on such factors as cloud top temperature, cloud 
growth, and departure of precipitable water from a summer- 
time normal. 

The Griffith/Woodley (Griffith et al, 1978) technique is 
designed to estimate rainfall in the tropics, over large 
space and time scales, using geosynchronous visual or 
infrared satellite imagery. This time-dependent technique 
was empirically derived as a relationship between cloud 
area, echo area, and rain rate for two areas in south 
florida, with raingauge-radar providing the ground truth, 
and was then tested in other tropical areas. This scheme 
was subsequently tested further in ex tr atropical areas 
(Griffith et al, 1980), with modifications to the rainfall 
amount predicted. 

The determination of a cloud area-rainfall relationship 
first required the specification of both a visual and an 
infrared threshold to define the cloud area. The visual 
brightness threshold, normalized for radiation geometry, was 
80 counts for the third Application Technology Satellite 
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(ATS-3) and the infrared threshold was 253 K (-20 C) . The 
thresholds were based on a comparison of the clouds with a 
given maximum digital count and the radar echoes associated 
with these clouds. 

The empirical cloud area-rainfall relationship was 
derived as a two step process. First, a relationship 
between the cloud area and the radar echo area, normalized 
for the maximum area achieved by the cloud or the cluster, 
was established for the visible and infrared satellite data. 
Second, the relationship between the echo area and the rain 
volume was determined and was of the form: 

Rv = lAe (A. 4) 

where fiv is the rain volume per hour (m3 -h-1), I is rain in 
units of (m3 -km-2 -h-1), and Ae is tne echo area (kn2) 
defined by the 1mm -h-1 rain rate. Thus, given a time 
sequence of convective clouds (or cluster areas) measured 
from visible or infrared satellite images, a volumetric rain 
rate can be estimated. 

Stout et al, (1979) modified the Griffith/Woodley tech- 
nique (Griffith et al, 1978) to estimate volumetric rain 
rate directly from a cumulonimbus cloud area and area change 
according to the equation: 

R = a(0)A + a(1)dA/dt (A. 5) 

where R is the volumetric rainfall of the cloud (m3 -s- 1) , A 
is the cloud area (m2) , dA/dt is the change of cloud area 
over time (m2 -s-1), and a (0) and a(1) are constants with 
dimensions (m -s-1) and (m) respectively. The two constants 
were calculated by a least square fit of cloud area-rain 
rate pairs based on visible and infrared geosynchronous 
satellite data and 5.3 cm ship radar rain data collected 
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during GATE. The cloud area and its change are defined by 
the threshold value. The visible threshold for cloud area 
calculations was 60 digital counts on the ATS-3 (corre- 
sponding to an albedo of 0. 45 with the sun overhead) , or 
eguivalently 172 digital counts on the Geosynchronous 
Heteorological Satellite (SMS 1) . The infrared threshold 
was 160 digital counts (-26 C) . The standard error between 
the estimated rainfall and the mean radar rainfall was 62% 
and 76% for the visual and infrared equations respectively. 

Wylie (1979) attempted to use the tropical convective 
rainfall techniques of Griffith et al, (1978) and Stout et 
al , (1979) for estimating precipitation in Montreal, Canada. 
Using visual satellite data, corrected for the changing sun 
angle (Mosher, 1975), infrared satellite data, and 10.0 cm 
radar measured rainfall rates, Wylie studied six days of 
precipitation, three days each in June and September 1977. 
Wylie concluded that because of air mass differences between 
Montreal and the tropics, the Griffith and Stout estimation 
techniques did poorly in Montreal, Canada. The singlemost 
important limitation with these two schemes was the diffi- 
culty in measuring cumulonimbus cloud area when the "anvils 
were often merged into large cloud masses and the extensive 
stratus cloud cover often obscured the pictures. " Wylie 
also noted that the Griffith et al, (1978) threshold of -26 
C had to be changed to -16 C for the summertime Montreal, 
Canada area. With the warmer cloud top temperatures the 
cloud area were a larger, more appropriate size for 
tracking. 

Wylie (1979) then attempted to combine sounding data 
input into a one-dimensional model (Simpson and Wiggert, 
1969) and satellite cloud cover measurements to estimate 
rainfall for Montreal. With the GATE measurements for rain 
rates associated with satellite- derived cloud areas and the 
model output, rainfall rates were estimated by multiplying 
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the two values. The most accurate estimations were for the 
cumulus clouds in the warm air masses occurring in June, the 
cases the model was designed to handle. Wylie concluded 
that in order to estimate the rainfall in all geographical 
areas and seasons a more sophisticated model would be 
needed. 

Negri and Adler (1981) did one case study of 15 thunder- 
storms in the Oklahoma, Arkansas, and ’lissouri area on 24 
April 1975. They used radar data for ground truth and had 
special 5 minute GOES-E satellite passes over the area of 
interest. They were able to determine that precipitation 
began falling, as indicated by radar data, for cloud top 
temperatures ranging from 229 K to 260 K (-44 C to -13 C) . 
The mean cloud top temperature value was 247 K (-26 C) . 



45 



appendix b 

aUENCH AND KEEGAN (1979) NORNALIZ ATION SCHEBE 



The Muench and Keegan (1979) normalization relates the 
normalized reflectivity (r^) to the varying solar angle and 
maximum digital counts through the reflectance term (r) , and 
the anisotropic scattering through the X term. Table XVII 
defines the symbols. Table XVIII lists the geometric iden- 
tity equations, and Table XIX lists the normalization equa- 
tions. Fig. 11 gives an example of the normalization 
applied to the stated location. 
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APPENDIX C 



FIGDBES 




Fig- 1. Geographical Location of Service-A Station 
Station Report Data 
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Fig. 2. No-Precipitation, Infrared Data Distribution 
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Fig. 3. No-Precipitation, Visible Data Distribution 
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Fig. 5- Precipitation, Infrared Data Distribution 
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Fig. 6. Precipitation, Visible Data Distribution 
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IR 




Fig. 7. Separation of Different Cloud Types in VIS 

and IB from TIEOS-N AVHRR-data. Sun Elevation 45 
(Liljas, 19 81) 




Fig. 8. Qualitative Indication of Precipitation 

Intensities Derived Using VIS and IH Data from 
TIEOS-N AVHRR-data. (Liljas, 1981) 
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Fig- 9. 



Fig. 10- 
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Frequency Plot of Rain Data Distribution for GATE 
day 248, 1300 GMT (Lovejoy and Austin, 1979) 
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Frequency Plot of No-Bain Data Distribution for 
GAT5 day 248, 1300 GHT (Lovejoy and Austin, 1979) 
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a 




Fig. 11. Normalized Cloud Reflectivity as a Function of 
Video Count and Time of Day for 16 April 1977 
at 42 N and 74 W with Calibration Count Equal 
to 60 (Muench and Keegan, 1979) 
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APPENDIX D 



Table I. 



TABLES 

Information Available for Each Observation 



1 . Station ID 

2. Latitude 

3. Longitude 

4. Altimeter 

5. Temperature 

6. Dew Point 

7. Wind Direction 

8. Wind Speed 

9. Sea Level Pressure 

10. Cloud Group 

11. Visibility 

12. Precipitation Amount 

13. Wind Gusts 

14. Current Weather 

15. Maximum Temperature 

16. Minimum Temperature 

17. Snow 

18. Observation Date/Time 

19. Visible/Infrared Satellite Data 
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Table II 



Compariscn of Means and Standard Deviations 
of Each Matrix Size tor Each Category 



No Precip 
II 


IQ. X 10 


8x8 


6x6 


4x4 


Mean 


275.82 


275.38 


275. 74 


275. 71 


St Dev 


18.37 


18.69 


19.03 


19. 37 


All Precip 

II 










Mean 


244.61 


244.00 


243. 47 


242. 99 


St Dev 


21.05 


21.26 


21.46 


21. 72 


No Precip 
VIS 










Mean 


26.91 


27.02 


27.15 


27. 29 


St Dev 


24.96 


25. 29 


25.72 


26.32 


All Precip 

II 










Mean 


53.69 


59.63 


60.5 7 


61.32 


St Dev 


27.15 


27.40 


27.90 


28. 53 
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Table III 



No-Precipitatioii/ Infrared Analysis 



Mean.' 275,71 St Dev: 19.37 



COLD TAIL: 1461 Obs ( < 240.0 K ) 





1 Qbs 


Hiah 


Middle 


Low 


Overcast • 


0.50 


0.35 


0. 10 


0.05 


Broken 


0.50 


0.4 1 


0.07 


0.03 


Total 


1,00 


0.76 


0. 16 


0.08 



Table IV. No-Precipitation, Visible Analysis 



Mean : 


27.29 


St 


Dev: 26.32 




BEIGHT 


TAIL: 999 


Obs 


( > 0.70 ) 






# Obs 


Hi^h 


Middle 


Low 


Overcast 


0,49 


0. 18 


0. 14 


0.17 


Broken 


0.51 


0.41 


0,07 


0. 03 


Total 


1,00 


0.59 


0. 22 


0.20 
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Table V. Precipitation, Infrared Analysis 



Kean; 242.99 St Dev: 21-72 





Cold 


.liddle 


Warm 




End 


Region 


Tail 


TEH 


0.15 


0. 05 


0.06 


TRW- 


0-42 


0.21 


0- 18 


RH 


0. 02 


0. 04 


0.02 


RW- 


0.20 


0.37 


0.50 


R 


0. 02 


0.01 


0.01 


R- 


0. 14 


0.29 


0.21 


TR 


0.00 


0. 01 


0.00 


TE- 


0.05 


0.03 


0.02 


1 Obs 


959 


626 


257 



Table VI. Precipitation, Visible Analysis 



Mean : 


6 1.32 


St Dev: 


28.53 




Bright 


Middle 


Dim 




End 


Region 


Tail 


TRW 


0- 18 


0. 1 1 


0.08 


TRW- 


0.29 


0. 29 


0.28 


RH 


0.06 


0. 03 


0.02 


EW- 


0. 19 


0.30 


0.43 


R 


0-03 


0.02 


0.01 


R- 


0. 21 


0. 22 


0.15 


TR 


0-01 


0. 01 


0.00 


TR- 


0.03 


0. 02 


0.03 


# Obs 


163 


916 


327 
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Table VII 



Comparison of Means and Standard Deviations for 
Precipitation with Broken/Overcast Cexlings 
with Paul (1983). (Paul’s values in 
parenthesis) 



Mean 
St Dev 



Precipitation 

11 X 13 

244. 6 K (252.8 K) 

21. 1 K ( 20. 5 K) 



- Infrared 

243.0 K (250.6 K) 
21.7 K ( 21.5 K) 



Mean 
St Dev 



Precipitation 
11 X 10 
0.587 (0.580) 

0.272 (0.210) 



Visible 

i X i 

0.613 (0.617) 

0.285 (0.225) 



Table VIII. Comparison of Means and Standard Deviations for 
No- Precipitation witn Broken/Overcast Ceilings 
with Paul (1983). (Paul's values in 
parenthesis) 



No-Precipitation - Infrared 

10 X 10 1 X 1 

Mean 275.8 K (279. 8 K) 275.7 K (279.7 K) 

St Dev 18.4 K ( 16.6 K) 19.4 K ( 17.8 K) 



Mean 
St Dev 



No- Precipitation 

10 X 10 

0.269 (0.271) 

0. 250 (0.189) 



- Visible 
4x4 

0.273 (0.277) 

0.263 ( 0.207) 
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Table IX 



Verification of SPADS Infrared Thresholds 



Obse rv ati ons: Yes 

Satellite 
Estimat e : 



Yes 

No 

Total Obs 



1195 (0.64) 
2505 (0. 13) 
3698 



No 



661 (0.36) 

16,35 1 (0. 37) 
17,512 



Total 

Obs 

1356 
19,354 
21,2 10 



Table X. Verification of SPADS Visible Thresholds 



Satellite 
m^e : 


Observations: Yes 


No 


Total 

Obs 


Yes 




909 (0.63) 


528 (0.37) 


1437 


No 




1962 (0.13) 


12,778 (0.37) 


14,740 


Total 


Obs 


2871 


13,306 


16,177 



Table XI. Verification of SPADS Visible and Infrared 
Thresholds Combined 



Obser vat ions: Yes 

Satellite 

Est im ate 



Yes 

No 

Total Obs 



343 (0.32) 
825 (0.07) 
1668 



No 



180 (0.18) 
10,319 (0.93) 
10,499 



Total 

Obs 

10 23 
11,144 
12,167 



Table XII 



Threshold Values Describing Precipitation 
Intensity Levels as Applied in Fig. 8 
(Liljas, 1981a) 



The Su m of Diaital Lev els 



Chi Ch 4 



291 - 310 Light Rain 

311 - 330 

331 - 350 

351 - 370 

371 - 390 

399 * Very Strong Rain 



Table XIII. Statistical Comparison of Rain Area Mapping 
Techniques (Lovejoy and Austin, 1979) 





Opt. 2-D Boundary 


IR Optimum Tlireshold 


Visible Optimum Threshold 


Arcs Rsin 


Total No. 
of Points 


Day 


(RmIR) * 100* 


1R(K) 


(Rm/R) X 100 


(Scale: 0-1) 


(RmJR) X 100 


Coverage (%) 


GATE 


242, 243. 246 
247.248, 231 
252. 261 


65 


<232 


53 


>0.68 


64 


15.8 


47706 


Montreal 


132 


56 


<232 


20 


>0.88 


31 


9.7 


40361 


180 


36 


<247 


55 


>0.80 


48 


24 0 


33738 


153 


53 


<254 


52 


>0.88 


49 


13.9 


22358 



•Referred to in text as ’’percentage of correct satellite rain.’* 



Table XIV. Statistical Comparison of the Accuracy of Rain 
Areas (Lovejoy and Austin, 1979) 



Technique 


Region 


Number of Images 
or Sequences 


Bias 


Error 

Factor 


^RMS 


2-D Pattern 
Matching 


Montreal 


17 


1.13 


1.26 


0.22 


2-D Pattern 
Matching 


Montreal 


3 


1.08 


1.19 


0.18 


Optimum IR 
Threshold 


Montreal 


3 


1.38 


1.74 


0.71 


Optimum Visible 
Threshold 


Montreal 


3 


1.54 


1.59 


0.58 


2-D Pattern 
Matching 


GATE 


8 


1.21 


1.41 


0.25 
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Table XV. General Coefficients for the Detemination of 
Precipitation Based on the Three Data Types 
(Wylie, 1982) 



CORE ELATION 

MEASURED WITH 6 HOUR NUMBER 



PAM METER PRECIP. REPORT 


OF CASES 


F 


Vertically integrated 
precipitafcle water vapor 


0.48 


196 


5 8* 


Cloud top brightness 


-0.44 


184 


44* 


Cloud top height 


-0.40 


190 


36* 


Moisture convergence 


0.38 


184 


3 1* 


Cloud top temperature 


-0.35 


199 


27* 


Bubble model predicted cond. 


0.27 


115 


9 


500 mb vorticity advection 


-0.21 


173 


8* 


Parcel lifted index 


-0.20 


200 


8* 


700 mb temperature advection 


0.20 


173 


7* 


Sfc temperature advection 


0.19 


156 


6 


850 mb temperature advection 


0.17 


189 


6 


Wind convergence (sfc) 


0.09 


167 


1 


Vertical wind shear 


0.03 


156 


0 



♦ Significant ccrrelation at the 99% level. 
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Table XVI. Summary of Life History Threshold Values 



STUDY 


LOCATION 


CASES 


TIME OF 
YEAR 


THRESHOLD 

INFRARED 


VALUE 

VISUAL 


Griffith, 
ct al. , 
(1978) 


Florida, 

Venezuela, 

Honduras, 

and hurricanes 

impacting East 

Coast United 

States 


34 days 


summers 

1969-1976 


-20®C 


80 counts* 


Stout, 
et al. , 

(1979) 


tropical North 
Atlantic 


57 obser- 
vations 


September 

1974 


-26°C 


0.45 albedo 
(sun over- 
head) 


Wylie 

(1979) 


Montreal 


6 days 


June 1977 
September 
1977 


-16®C 


- 


Negri and 
Adler 
(1981) 


Oklahoma , 
Arkansas, 
Missouri 


1 day 

(15 thunder- 
storms) 


April 24, 
- 1975 


-27°C 


- 



* ATS-3 satellite 
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Table XVII. List of Symbols 



Symbol 


Description 


Units 


c 


GOES video count number (0-63) 


dimensionless 




GOES video count number for perfect 
diffuse reflector and overhead sun 


dimensionless 


G 


Greenwich meridian time 


hours-minutes- 

seconds 


R 


Distance of earth to sun 


km 


% 


Mean distance of earth to sun 


km 


d 


Julian date 


dimensionless 


Ii 


Hour angle 


radians 


r 


Cloud reflectivity 


dimensionless 


Y 


Arc-length observer to subsatellite 
point 


radians 


S 


Declination of the sun 


radians 


C 


Zenith angle of the sun 


radians 


A 


Longitude 


radians 


As 


Longitude of subsatellite point 


radians 


X 


Anisotropic scattering coefficient 


dimensionless 


♦ 


Latitude 


radians 


h 


Azimuth of the sun 


radians 


^2 • 


Azimuth of the satellite 


radians 
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Table XVIII 



Basic Geometric Satellite- Earth Relationships 



* 



Declioation : 

6 - 0.408 sin [(d-81) • 2ir/365] 

Solax distance ratio: 

R/Rq " 1 “ 0.167 cos I(d-14) * 27r/365] 
Hour angle: 

h » A IT - G(hours) x it/ 12 
Arc-length: 

COSY ■ cosCAg - A) cos4 
Satellite azimuth: 

sin C <>2 - If) ■ sin (Ag - A)/sinY 
Solar azimuth angle: 

cose ■ sin^ sin6 + cos4 cos6 cosh 
Solar azimuth: 

sin^j^ - cos<S sinh/sine 



Angles in radians 

Table XIX. Muench and Keegan (1979) Normalization 
Equations * 

f = • sec? 

” I’f’2 ” 

= cos^ CC? - 50) * 1 . 8 ) 

= 0.7 cos C(C - 22.5) * 4) ♦ (1 - cos? ) 

C 3 = cos® CCA(() - 70) » 1,3) 

X = 1.0 + 0.05 * (1 + cos(2»?)) + 0.20 (C^ + C 2 ) C 3 

“ U.09 - 2*(1.09 - r * X *(R/Rq)^)/( 1 + cos^'^^ c)] 

★ 

Angles in degrees 
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